Characterizing the Peak in the 
Cosmic Microwave Background Angular Power Spectrum 
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A peak has been unambiguously detected in the cosmic microwave background (CMB) angular 
spectrum. Here we characterize its properties with fits to phenomenological models. We find that 
the TOCO and BOOM/NA data determine the peak location to be in the range 175-243 and 151- 
259 respectively (both ranges 95% confidence) and determine the peak amplitude to be between 
« 70 and 90 /j,K. By combining all the data, we constrain the full- width at half-maximum to be 
between 180 and 250 at 95% confidence. Such a peak shape is consistent with inflation-inspired flat, 
cold dark matter plus cosmological constant models of structure formation with adiabatic, nearly 
scale-invariant initial conditions. It is inconsistent with open and defect models. 
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Introduction. If the adiabatic cold dark matter (CDM) 
models with scale-invariant initial conditions describe our 
cosmogony, then an analysis of the anisotropy in the 
CMB can reveal the cosmological parameters to unprece- 
dented accuracy [Q . A number of studies have aimed at 
determining, with various prior assumptions, a subset of 
the ~ 10 free parameters that affect the statistical prop- 
erties of the CMB |^||]. The parameter most robustly 
determined from current data is Q, the ratio of the mean 
matter/energy density to the critical density (that for 
which the mean spatial curvature is zero). These in- 
vestigations show that f2 is close to one. This result, 
combined with other cosmological data, implies the ex- 
istence of some smoothly distributed energy component 
with negative pressure such as a cosmological constant. 

A weakness of previous approaches |^|| is that the con- 
clusions depend on the validity of the assumed model. In 
this Letter we take a different tack and ask what we know 
independent of the details of the cosmological model. We 
find the peak location, amplitude and width are con- 
sistent with those expected in adiabatic CDM models. 
Furthermore, as ^ pca k — 200S1 -1 / 2 in these models, the 
observed peak location implies Q ~ 1. The determina- 
tion of the peak location is robust; it does not depend on 
the parametrization of the spectrum, assumptions about 
the distribution of the power spectrum measurement er- 
rors, nor on the validity of any one data set. The model- 
dependent determinations of f2 are further supported by 
the inconsistency of the data with competing models, 
such as topological defects, open models with < 0.4, 
or the simplest isocurvature models. 
The Data. The last year of the 1000's saw new results 
from MS AM f§, PythonV g, MAT/TOCO §,§, Viper 
@, CAT §, IAC @ and BOOM/NA 0, all of which 
have bearing on the properties of the peak. These results 
are plotted in Fig. 1. We have known for several years 
that there is a rise toward towards I = 200 but it is now 



clear that the spectrum also falls significantly towards 
I = 400. 

For all the medium angular scale experiments, the 
largest systematic effect is the calibration error which 
is roughly 10% for each. Contamination from foreground 
emission is also important and not yet fully accounted 
for in some experiments (e.g. TOCO). A correction for 
this contribution, for which STi ~ l^ 1 / 2 , will affect the 
amplitude of the peak though will not strongly affect 
its position. Thorough analyses by the MSAM |L2| and 
PYTHON || teams show that the level of contamination 
in those experiments was < 3%. 

The three experiments that have taken data that span 
the peak are MSAM, TOCO, and BOOM/NA. All exper- 
iments exhibit a definite increase over the Sachs- Wolfe 
plateau though the significance of a feature based on the 
data alone, e.g. a peak, differs between experiments. We 
may assess the detection of a feature by examining the 
deviation from the best fit flat line, ST. For the three 
MSAM points, we find ST = 46 ± 4.9 with a re- 
duced x 2 °f 0-43 (Probability to exceed, P >x z = 0.65. 
The calibration error is not included.). Thus, no fea- 
ture is detected with these data alone though there is a 
clear increase over DMR (l|. For the seven BOOM/NA 
points, we find ST — 55.3 ± 4.2 /xK with a reduced 
X 2 of 1.94 (-P> x 2 = 0.05, assuming the data are anti- 
correlate^ at the 0.1 level Q). For the ten TOCO 
points, ST = 69.3 ± 2.7 with a reduced x 2 of 4.86 
{P> x 2 < 10~ 5 ) Calibration errors will not change x 2 l v i 
however a correction for foreground emission will have 
a slight effect. Though we examine all data in the fol- 
lowing, we focus particularly on BOOM/NA and TOCO 
because of their detections of a feature. 
Fits to Phenomenological Models. To characterize the 
peak amplitude and location we fit the parameters of two 
different phenomenological models. For the first, we start 
with the best fit DK99 @ adiabatic CDM model, STf 
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and form STi = (ST l DK - ST^ Q )a + ST^f Q by varying 
a, and then stretching in I jl3]. We characterize each 
stretching with the peak position and peak amplitude. 
This method has the virtue that the resulting spectra 
resemble adiabatic models and so if one assumes that 
these models describe Nature, then these results are the 
ones to which we should pay the most attention. 
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FIG. 1. Bandpowers from TOC097 (cyan open trian- 
gles), TOC098 (blue filled triangles), BOOM/NA (green 
filled squares), MS AM (red open squares), CAT (black 
open pentagon), I AC (black filled pentagon), PyV (black 
open circles) and Viper (green filled circles). The y-axis 
is STi = + l)Gi/(2w) where Ci is the angular spec- 

trum. The models are, peaking at left to right, the best 
fit models of M for Q, = 1, fi = 0.4 and Q = 0.2. The 
Q, = 1 model has a mean density of non-relativistic matter, 
fim = 0.31, a cosmological constant density of J1a = 0.69, 
a baryon density of fit = 0.019ft -2 jl5), a Hubble constant 
of Ho = 65 km/sec/Mpc, an optical depth to reionization of 
r = 0.17 and a power spectrum power-law index of n — 1.12, 
where n = 1 is scale-invariant. The shaded areas are the re- 
sults of fitting the power in 14 bands of I to all the data (from 
1999 and previous years) as in Jli| . Many of the bands are at 
low I and cannot be discerned on this plot. Calibration errors 
are not shown though are included in the best fit. 



Our second model for 



SJf 



is a 



ST? = 



Gaussian 

A 2 exp ^— (I — l c ) 2 /(2of)j. Depending on the width, 
this spectrum can look very much like, or unlike, the 
spectra of adiabatic models |L7| . We view this versatility 
as a virtue since we are interested in a characterization 
of the peak which is independent of physical models. 

We fit to these phenomenological models in two ways. 
For the stretch model, we examine the \ 2 °f the residuals 
between the published data and each model. The widths 



of the window functions are ignored and we assume the 
data are normally distributed in 5Ti with a dispersion 
given by the average of the published error bars (GT 
in Table 1). This is an admittedly crude method but it 
works well because the likelihoods as a function of STi are 
moderately well approximated by a normal distribution. 

For both the Gaussian shape and the stretch model, 
we also perform the full fit as outlined in BJK |ll| (RAD 
in Table 1). For the Gaussian shape model, the con- 
straints on the amplitude and location are given below 
after marginalization over the width <r;. In all fitting, we 
ignore the experiments that are affected by I < 30 (DMR, 
FIRS [|l8| and Tenerife 19 1) because we want the param- 
eters of our Gaussian to be determined by behavior in 
the peak region. 
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0.85 


206 ± 26 
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33/30 1.13 


0.28 


262 ± 24 


68 



*ALL stands for all publically available data sets (except 
for VIPER which was not used because of unspecified 
point-to-point correlations), the T is for the TOCO data, 
the B for BOOM/NA and the P is for "Previous" , mean- 
ing all data prior to BOOM/NA and TOCO. f G and S 
are for the Gaussian shape and stretch methods respec- 
tively *iV is the number of data points and v the degrees 
of freedom. §Rado and Radoo corresponds to log normal 
and normal distributions for the likelihood respectively. 

The main thing to notice in the Table is that the posi- 
tion of the peak is robustly determined by either TOCO 
or BOOM/NA to be in the range 185 to 235, regardless 
of the method. For the quoted errors, we have marginal- 
ized over all parameters except the position. The peak 
amplitudes are subject to change as there is some depen- 
dence on the model parametrization and the foreground 
contamination has not been thoroughly assessed. 

We account for the calibration uncertainty through a 
convolution of the likelihood of the fits with a normal 
distribution of the fractional error (HQ. BOOM/NA, 
TOC097 and TOC098 have calibration uncertainties 
of 8%, 10.5% and 8% respectively. However, 5% of 
this is due to uncertainty in the temperature of Jupiter 
and therefore, assuming that these uncertainties add in 
quadrature, we get erj up = 0.05, <JT97 — 0.092, (7T98 = 
0.062 and <j B 97 = 0.062. We then find, for TOCO, that 
the full likelihood in STi and I is given by 
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L(l c ,5Tl) = J d<Tiduj up du T 97duT98L T 97(lc,STiU,i up UT97,Vl) 

xL T9 $(l c , 5TiUj up u Tg8 , ai)P G (u 97 - 1; 07-97) 

xP G (u Tg8 - 1; CTT98 )Pg{u Jup - 1; 0J U p) (1) 

where Pc(x;a) = exp (— x 2 /(2er 2 )) /\/27rcr 2 , u is in- 
tegrated from to oo and, e.g., LT97(l c ,STi,ai) — 
cxp(— x 2 /2) where x 2 is evaluated on a grid of <5T 2 , l c 
k 01 using RADPACK [fHJ as discussed in BJK. We get 
similar results for TOCO when simply using a combined 
total calibration error of 8.5%. 

For the Gaussian model we can also marginalize over 
A and l c to place 95% confidence bounds on the width: 
75 < a t < 105 for ALL, 50 < a t < 105 for TOCO and 
55 < ai < 145 for BOOM/NA. 
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FIG. 2. Likelihood contours for I vs STi for the position 
of the peak. For BOOM/NA and TOCO, we use the stretch 
method using RADPACK J2l| and include the calibration er- 
ror. For Previous and ALL (tightest contours) we have not 
used generalizations of Eq. [l], but instead have fixed the cal- 
ibration parameters. All contour levels correspond to 5%, 
68%, and 95% enclosed, or roughly the peak, la, 2a. 

Are the data in Fig 1 consistent? DK99 found that the 
best-fit model, given all the data at the time, had a \ 2 of 
79 for 63 degrees of freedom, which is exceeded 8% of the 
time. Here we see that the % 2 for the fit of the Gaussian 
model is 69 for 55 degrees of freedom, which is exceeded 
10% of the time. We conclude that, although there may 
well be systematic error in some of these data sets, we 
have no compelling evidence of it. However, we take cau- 
tion from the fact that we had to adjust the calibration 
parameters from their nominal values to their best-fit val- 
ues in order to reduce the \ 2 to 69. Left at their nominal 
values with calibration uncertainty ignored, the data are 
not consistent with each other. Thus we believe that the 



compilation results are perhaps less reliable than those 
for either BOOM/NA or TOCO. 

Implications for Physical Models. Flat, adiabatic, nearly 
scale-invariant models have similar peak properties to 
those of our best-fit phenomenological models. Most im- 
portantly the peak location, as determined by three in- 
dependent data sets ("Previous", TOCO, BOOM/NA), 
is near I ~ 210, as expected. Depending on the data set 
chosen, the amplitude is higher than expected but can 
easily be accommodated, within the uncertainties, with 
a cosmological constant. Combining all the data, there 
is a preference for Z poa k > 210 which suggests a cosmo- 
logical constant |22j (at h — 0.65, / poa k goes from 200 at 
Qj\ = to 220 at £1a = 0.7). However, this result is not 
seen in any individual data set. 

A good approximation to the first peak in the DK99 
best-fit model is given by the Gaussian model with 07 = 
95. From the 07 constraints quoted earlier we see that 
the data have no significant preference for peaks that 
are either narrower or broader than those in inflation- 
inspired CDM models. 

A general perturbation is a combination of adiabatic 
and isocurvature perturbations. Adiabatic perturbations 
are such that at each point in space, the fractional fluc- 
tuations in the number density of each particle species is 
the same for all species. Isocurvature perturbations are 
initially arranged so that, despite fluctuations in individ- 
ual species, the total energy density fluctuation is zero. 
Given multiple components, there are a number of differ- 
ent ways of maintaining the isocurvature condition. Be- 
low we assume the isocurvature condition is maintained 
by the dark matter compensating everything else. 

Isocurvature initial conditions result in shifts to the 
CMB power spectrum peak locations. For a given 
wavenumber, the temporal phase of oscillations in the 
baryon-photon fluid depends on the initial relation be- 
tween the dark matter and the fluid. Those waves with 
oscillation frequencies such that they hit an extrcmum 
at the time of last-scattering in the adiabatic case, will 
hit a null in the isocurvature case (2^] . The effect on the 
first peak is a shift from I ~ 200ft" 1 / 2 to I ~ 350ft" 1 / 2 . 
Given the observation of ^ poa k — 210, simple isocurva- 
ture models require ft > 2 — which is inconsistent with a 
number of observations EG] . 

Critical to the Doppler peak structure, in either adia- 
batic or isocurvature models, is the temporal phase co- 
herence for Fourier modes of a given wavenumber [p4[ . 
In topological defect models, the continual generation of 
new perturbations by the non-linear evolution of the de- 
fect network destroys this temporal phase coherence and 
the acoustic peaks blend into a broad hump which is 
wider and peaks at higher I than the observed feature. 

One can make defect model power spectra with less 
power at I = 400 than at I = 200 with ad-hoc modifica- 
tions to the standard ionization history [ B5j . But even for 
these models the drop is probably not fast enough p7|. 
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The contrast between the power at I — 200 and I = 400 

is a great challenge for these models. 

T here are scenarios with initially isocurvature condi- 
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tions that can produce CMB power spectra that look 
much like those in the adiabatic case. This can be done 
by adding to the adiabatic fluctuations (of photons, neu- 
trinos, baryons and cold dark mater) another compo- 
nent, with a non-trivial stress history, which maintains 
the isocurvature condition p^j. 

Conclusions. Our phenomenological models have allowed 
for rapid, model-independent, investigation of the con- 
sistency of CMB datasets, and of the robustness of the 
properties of the peak in the CMB power spectrum. The 
peak has been observed by two different instruments, and 
can be inferred from an independent compilation of other 
data sets. The properties of this peak are consistent with 
those of the first peak in the inflation-inspired adiabatic 
CDM models, and inconsistent with a number of com- 
peting models, with the possible exception of the more 
complicated isocurvature models mentioned above. It is 
perhaps instructive that where the confrontation between 
theory and observation can be done with a minimum of 
theoretical uncertainty, the adiabatic CDM models have 
been highly successful. 
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